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One of the most promising ways of detecting primordial gravitational waves generated during 
inflation is to observe B-modes of polarization, generated by Thomson scattering after reionization, 
in the cosmic microwave background (CMB). Large scale foregrounds though are expected to be a 
major systematic issue, so - in the event of a tentative detection - an independent confirmation of 
large scale gravitational waves would be most welcome. Previous authors have suggested searching 
for the analogous mode of cosmic shear in weak lensing surveys but have shown that the signal to 
' noise of this mode is marginal at best. This argument is reconsidered here, accounting for the cross- 

ed ' correlations of the polarization and lensing B-modes. A lensing survey can potentially strengthen the 

(— I , argument for a detection of primordial gravitational waves, although it is unlikely to help constrain 

' the amplitude of the signal. 

. PACS numbers: 95.35.+d; 95.85.Pw 

(N 

o, 

. Many models of inflation predict that gravitational waves were produced along with density perturbations in 
(-H ' the early universe [ll-IE|- While the density perturbations have been detected and mapped out very precisely, the 
O ^1 gravitational waves have yet to be detected. Perhaps our best hope is to observe the unique pattern of polarization 
' ■ - so-called B-modes - that gravitational waves produce when photons in a warped radiation field scatter off free 
electrons ^-Q- One exciting feature of this signal is that it was produced at two distinct epochs: first, during 
recombination when the scattering rate starts to tail off and photons can travel longer distances, and, then, after 
^ I reionization at much later times. These two epochs are imprinted in the angular spectrum on large scales (multipoles 
'— '. 1^2 — 10) and small scales (~ 50 — 100). Hence primordial gravitational waves leave the distinct signature of a 

^ ' double-humped peak in the power spectrum of polarization B-modes. 
^ Detection, however, is far from guaranteed. Among the hurdles that must be overcome are: (i) unknown amplitude 

of the signal [Ij, (ii) high sensitivity [11] needed to detect even the largest signal possible, and (iii) foreground 
7-H , contamination |12|. Even if the first two of these are overcome, and a B-mode detection is claimed, distinguishing a 
putative signal from foregrounds is likely to remain a significant problem. Here, I explore one possible way to test a 
future B-mode detection against foreground contamination: supplementing the polarization B-mode signal with an 
^-H . all-sky, deep weak lensing signal. 

' The idea is simple: gravitational waves also leave an imprint on the cosmic distortion tensor that governs the 
, propagation of light over long distances, and hence the observed shapes of galaxies. Density perturbations in the 
' distortion tensor produce scalar shear, which was first detected in 2000 |13l4l6| and has since been measured over 
J> , larger scales by many experiments (see, e.g., [l3| for a review), but gravitational waves produce a different type of 
distortion [l^, sometimes called the pseudo-scalar shear, or the curl mode of shear, the exact analogue of the B-mode 
of the polarization field. A number of studies [Tol . [20| have shown that the power spectrum of this mode, even in 
■ the most optimistic scenarios (large amplitude, all-sky survey with low shape noise), will be difficult to detect except 
" " ' perhaps at the lowest multipoles. Here I point out that the situation is not quite as grim as these studies suggest. 
First, even a moderate signal to noise independent detection in an arena with completely different systematics would 
significantly increase our confidence in a putative CMB detection. Further, the cross-correlation between the CMB 
and lensing B-modes is non-zero and this helps the search for a detection. 

The main purpose of this paper is to calculate the extent to which this is true, the extent to which the same modes 
which produce the CMB B-mode reionization signal are responsible for the B-mode in a lensing experiment. This 
is an interesting academic question in its own right: two very different phenomena owe their existence to the same 
spectrum of gravitational waves, and it is interesting to quantify how correlated they are. Sections H and HI present 
the formula for the auto-correlations of the B-modes in these two different arenas. Since I am interested only in the 
large scale feature, the focus in §H is on the reionization signal, and, as far as I can tell, the final formula for this signal 
in Eqs. and (IT5|) has not been presented elsewhere. In §IV, I calculate the cross-correlation between the two and 
show that it is non- negligible. Finally, §V addresses the more general question of how much additional information 
can be gleaned from a lensing survey, taking into account the cross-correlation. Some of the details of the calculations 
are relegated to appendices. 
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There is one scenario for which the probe introduced here would be of tremendous utihty. Over a decade ago, 
cosmologists studied gravitational wave production in open inflation models [2l| - |25| and concluded that the amplitude 
would be boosted near the horizon. This would lead to a much different spectrum than the double humped peak which 
modulates the scale-invariant standard inflation spectrum. Recently, there have been suggestions [1^ [23| that inflation 
was preceeded by an earlier epoch of vacuum domination such that our Universe was produced via a Colcman-De 
Luccia tunneling process (28| and is therefore open, albeit with a small value of the curvature density. Even with the 
small curvature expected in our Universe, though, the large scale gravitational wave spectrum might be different than 
the simple scale-invariant one generally assumed. A test - such as the B-mode of cosmic shear - sensitive to only the 
largest moments seems tailor-made to probe this class of models. 



II. B-MODES OF POLARIZATION 



The CMB radiation is described by a 2 x 2 intensity tensor (see, e.g., Ref. [2^, |30|), the traceless, symmetric part 
of which consists of two fields Q{h) and U{h). It is useful to decompose Q and U into E/B modes, where the B-mode 
multipole moments are given by 
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<fn [(2l^:„(n)) (g + iU){n) - (_2l^:„(n)) {Q - iU){n)\ 



(1) 



where the ±2^Zm are the spin-2 spherical harmonics. The polarization fields Q{n) and U{'h) are induced by Thomson 
scattering in the presence of a quadrupole distribution, so the polarization produced after reionization is 



(Q ± iU) (h) = 



5^6 Jo 



-Ore 



(2) 



m" = -2 



where 77 is the conformal time associated with the comoving distance D (in the flat universe we will assume throughout, 
T] = rio — D). That is, one integrates along the line of sight out to a spherical shell a distance Droion away (see Fig. 
[U adapted from Ref. [3l|), weighting by the scattering rate f = ngaTa (where Ue is the free electron number density, 
tJT the Thomson cross-section, and a the scale factor), and the different components of the quadrupole 92m- Eq. ([2]) 
relates the polarization field produced after reionization to the quadrupole of the radiation field, no matter what the 
origin of that quadrupole. The quadrupole produced by gravitational waves (GW) is of course linearly related to the 
tensor perturbation hij{x,r]) with (e.g., Ref. |32|) 



AD 



d(AD') 



—^n^n^hij{x + x; rj') 



(3) 



where AD = — is the distance between the point of interest and the surface of last scattering, h = dh/drj', and 
x' = AD'fi' identifies a point within a sphere of radius AD centered at x (see Fig. [Ij. To be clear, the argument 
77' in the tensor perturbation denotes the conformal time when the GW (which travels at the speed of light) was at 
position X + x' . That is, it was a distance AD' away from the point [x = (Dh)]. This corresponds to conformal time 
ri' = r]- AD' = rjo - (D + AD'). 

The tensor field is the sum of two independent modes in Fourier space: 



hij{x;r]) 



(2^ 



a— 



(4) 



where T(fc, rj) captures the evolution of the GW when they enter the horizon. In a matter dominated universe, 
T{k,ri) — 3ji{kr])/{kr]), but we use the exact solution to the transfer function here which differs slightly at late times 
due to dark energy domination. The orientation of the two modes depends on the direction of the k vector. If k is 
chosen to lie along the z-axis, then 



.(+) 





(5) 



More generally, they are transverse, traceless matrices normalized so that Tr[e(")e^'^)] = 2Sa 
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FIG. 1: The double integral which determines the B-mode signal due to reionization. The outer integral (Eq. ((2])) is along 
the line of sight Dh out to the time/distance of reionization, Droion- The quadrupole in the integrand is determined by an 
integral along AD'h' (Eq. ([S])) out to the surface of last scattering (a distance AD away) to capture the contribution of the 
gravitational waves. 



Armed with these results, we can write 



^-^/tSs E h^''\k)TrJ''\k) (6) 



where the transfer function for the B-modes in polarization is 



T^L^'^^^) - -'\hi^J J2 / d{AD')T{k,rj')jl:hk,ADy, (7) 



and 

Iimrn"ik,D) EE J d^ne*"^ [ (2^^ (n)) (2^2™" (")) - (-2l^:„(n))(-2r2,„.(n))] . (9) 

The transfer function is difficult to compute for general k but simplifies considerably when k lies along the z-axis. 
Fortunately, the auto- and cross-spectra that can be observed are rotationally invariant so in the end we will need 
only this simple case. The calculation is presented in Appendix |A] with the result that 



^"d(A7^')T(fc,,')f^. (10) 
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The m = 2 component is plotted for the lowest moments in Fig. [2] for Zioion = 10. 
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FIG. 2: The m = 2 moment of the B-mode transfer function in polarization for Z = 2 — 6 as a function of wavenumber k. Here 
reionization is assumed to occur instantaneously a.t z = 10. 

We will use Eq. (fTO|) to compute the cross-spectra with lensing. But, as long as we have it, we can first use it to derive 
a simple formula for the auto-spectrum of polarization B-modes due to reionization. Of course this spectrum is by now 
a standard feature of freely available codes which compute temperature and polarization two-point functions [s^, H^l , 
but it is nice to capture the physics in a simple semi-analytic formula. The amplitude of the GW mode is drawn from 
a Gaussian distribution with the same power spectrum: 



In standard slow roll inflation, the power spectrum is nearly scale invariant: Ph{k) / {2tt'^) ~ (4/7r)(i7//r7ipi)^, where 
Hi is the expansion rate during inflation. Squaring Eq. ^ and taking the expectation value using Eq. then leads 
to the polarization B-mode spectrum 



{U-^^kYh^fi^^k')) = {2^)H^05\k-k')PH{k). 



(11) 




m— — l 
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where 
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X [{l + 2)ji_^[kD)-{l-l)]i+^{kD) 



d{AD')t{k,r]') 



j2{kAD') 
{kAD'y 



(13) 



III. B-MODE OF COSMIC SHEAR 



The deformation tensor which describes propagation of light through the inhomogeneous universe can also be written 
as a 2 X 2 matrix, and its traceless, symmetric part can also be characterized by two fields: the two components of 
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shear 71 and 72 replacing Q and U . Thus, the moments of the B-mode of cosmic shear can be defined just as in Eq. ([TJ. 
There is a difference between the full deformation matrix and the full intensity matrix. In addition to Q and [/, the 
polarization matrix contains a piece describing the intensity / and circular polarization V . Neither / nor V is related 
to Q and U (in cosmology, / of course describes the CMB temperature anisotropics while V is thought to vanish 
because Thomson scattering induces no circular polarization). The analogous quantities in the lensing deformation 
tensor are the convergence k, and the rotation w. These are related to the two components of the shear, with k equal 
to the i?-mode and to to the i?-mode. Stebbins fis'l first examined the structure of this matrix and derived a number 
of useful relations among its components. Starting from the equivalent of Eq. ([T|) (with 71 and 72 replacing Q and 
U), one arrives at a simple expression for the i?-mode moments in terms of the rotation field 

Bi„r = j <fnYC^{n)uj{h). (14) 

The rotation receives no contributions from scalar perturbations, so it is non-zero only if tensor modes are present 
(at first order; see Ref. [131 for second order scalar effects). Explicitly, the moments of the B-mode of the shear field 
due to GW's are [H 

^"»= H h^-\k)Tt'-\k) (15) 

with the lensing transfer function defined as 

Tt^"\^)^\ j d^nY;^{n)e,,,ffh'e^^\k)k, J^De^^'^^^^ (16) 

where here Cijk is the 3D Levi-Civita symbol as opposed to the polarization tensor e''^\k). The integral here is out 
to the source galaxies, all assumed to be at distance Dg- 

Again the integrals over angles simplify considerably when k is chosen to lie along the z-axis. Appendix IB] contains 
the details leading to 

rV+lTA A 1 AG + 2)(/ + 1)^(^-1) ,^ .r^rr^n. n 7^ 

Ti,^' '{kz)^{i)+ [5m,2-5,n,-2\^^l + l ^ k dDT{k,r]o- D)jj-j^. (17) 

Fig. |3] shows the transfer functions for the lowest multipoles, assuming all galaxies are at redshift Zs = 1. The 
amplitudes here are important: the transfer function for Z = 2 is largest, reflecting the fact that the signal will be 
largest on the largest scales. 



IV. CROSS-CORRELATION 



We can now collect the window functions for the two auto-spectra and the cross-spectrum and compute the corre- 
lation co-efficient. Explicitly, in addition to Eq. (fT3|) . we have 



Wl^'^ik) = Tr{2l + l){l + 2){l + l)l{l-l)) 



dDT{k, fjo - D 



jiikD) 



[kDf 



Wf^{k) = &'K^{l + 2){l + l)l{l~l)k / dD' T{k, T]o - D') 



,,MkD' 



XT (ryo - D) [(/ + 2)ji^,{kD) -{I- l)ji+iikD)] 



{kD'f 

AD 



dD 



d(Ai^')T(fc,,y^('^^'^ 



(fcAL>')2 

The correlation coefficients which express the degree to which the moments are correlated are defined as 



(jPPfjLL 



(18) 



(19) 



where each Ci is an integral over the power spectrum modulated by the window function Wj. The window functions 
for the lowest moments of the cross-spectra are plotted in Fig. 2] when the background galaxies in the lensing survey 
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k/Ho 

FIG. 3: The m = 2 moment of the B-mode lensing transfer function for ? = 2 — 6. Sources are all assumed to be at z = 1. 




k/H„ 

FIG. 4: Cross-correlation window function for lowest moments when lensing galaxies are all aX z = 1 and reionization occurs 
instantaneously at ^rcion = 10. 

are all assumed to be at z = 1 and reionization takes place instantaneously at z — 1. As is clear, the I — 2 
moment has a negative correlation coefficient and is particularly sensitive to modes of precisely the size of the 
horizon. The lowest correlation coefficients - again when Zsouicc = 1 and Zioion = 10 - are: (02,03,04,05,00 = 
-0.32, 0.10, 0.31, -0.09, -0.20). 

A fascinating aspect of the cross-correlation is that, because of the oscillations in the gravitational waves, the 
correlation coefficient is sensitive to the redshift of the source galaxies in the lensing survey. Fig. [5] shows this 
dependence. In principle, one could imagine exploiting this dependence by weighting each galaxy by some factor to 
maximize the signal to noise extraction, but I do not pursue this possibility here. 



7 




V. SIGNIFICANCE 



How much would a lensing survey sensitive to the lowest moments add to our knowledge about the primordial 
gravitational waves? And how does the cross-correlation affect this answer? As political pollsters have discovered, 
the answers to these questions depend sensitively on how they are worded. 

For example, a simple question one might ask is how much tighter the constraints on the gravitational wave 
amplitude would become if the information from an all-sky lensing survey were added to that from a CMB polarization 
experiment. A quantitative answer to this can be obtained with the Fisher formalism. Consider two data points, the 
I = 2,m = 2 moment of the CMB B-mode and the same moment of the lensing B-mode. Normalize each so that the 
expected contributed to the variance from noise is unity. Also, allow for one free parameter, the normalization of the 
GW amplitude. A, with "true" value equal to 1. Finally, call the ratio of the signal variance to the noise variance A 
for each probe. Then the 2x2 covariance matrix is 



„_/ AX p + 1 aAy/XpXL \ 
^-[aA^X^ AXl + 1 )■ 

Here a is the correlation coefficient computed in the previous section. The n entry contains contributions from both 
signal {aXp) and noise (1). The signal contribution to the variance is Xp = CCJ^/Nj^ where C^^ is the power 
computed in §H and N^^ is the noise variance which of course depends on the experiment. Here and throughout this 
section, I will focus only on the lowest I — 2 moment since lensing contributions fall off rapidly with I. For orientation, 
current bounds on the gravitational wave amplitude from the CMB restrict CfJ^ to be of order (100 nK)^ while future 
experiments [s^ aim for noise as low as N^^ ~ (InK)^, so thoughts of Ap as large as ~ lO** are not crazy. Our goal 
is to estimate the uncertainty on the one free parameter A. The off-diagonal elements are free from noise so depend 
only on the cross-correlation of the two signals. Starting from this covariance matrix, and including all five I — 2 
moments, the standard formula for the Fisher matrix from a Gaussian process leads to a fractional error on A of 

_ P2 l + AL + Ap + ApAL(l-a') 

A V 5 [2A|,A2 (1 - a2)2 + 2(1 - a2)ApAi(Ap + A^) + A^ + A|, + 2a2ApAL]'/' ' 

The factor of -s/S in the denominator here comes from summing all five I = 2 harmonics. 

In the limit that Ap is large (corresponding to high signal to noise detection in the CMB), Fig. |6] shows how the 
constraints on the amplitude depend on the correlation coefficient and the lensing signal to noise. When the lensing 
signal is very small (A^ = 0.1) the error is simply equal to -^2/5 = 0.63, the minimum possible error due to the cosmic 
variance of the I = 2 CMB mode. If the signal to noise of lensing is higher, then the fractional error can go down. 
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in principle as low as a/I/S since the amount of information doubles. Fig. [6] illustrates that correlations degrade the 
extraction of the amplitude. This makes sense: when trying to measure a variance one wants as many independent 
numbers as possible. If a = 1, then the two numbers (B-modes from CMB polarization and from lensing) are not 
independent so the error on the variance increases. As a — >■ 1, the constraint on A reverts back to the \/2/5 limit 
that would be obtained without the additional lensing information. 

The conclusion from this exercise is that lensing information would help very little in the effort to pin down the 
gravitational wave amplitude. At best - information only from / = 2 and large signal to noise from both sets of 
experiments ~ the reduction in the error would be only a factor of -\/2- Most likely, if GW are detected, the I > 2 
moments in the CMB will be very important while those higher moments will be undetectable in lensing. So the 
gain from a lensing survey would be diluted significantly. The effect studied here - cross-correlations between the 
two signals - serves to further dilute the impact of lensing, as the information would be redundant and hence useless 
in constraining the GW amplitude. However, the results of the previous section suggest that is likely to be small 
enough so that the dilution would be minimal. If the lensing signal were detected, it would provide - for the most 
part - independent information about the GW amplitude. 
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FIG. 6: Constraints on the gravitational wave amplitude as a function of the cross-correlation co-efEcient if the signal in the 
CMB experiment is large. The signal to noise in the lensing survey is Xl ■ For large values of Xl , the error on the gravitational 
wave amplitude goes down by \/2 unless correlations between the two sets of measurements are large (a — >■ 1). 



Another way to probe the importance of lensing and cross-correlations is to focus not on parameter determination 
(A^) but rather on the firming up the case for detection. To show the difference between these two sets of ques- 
tions, consider a simple example: temperature anisotropies in the I = 2 mode as measured by WMAP [s^. The 
expected value of l{l + l)C//(27r) = 3C2/7r in the standard ACDM model is 928{iJ,K)'^, while the measurement error 
contributes only 0.01{fj,K)^ , so the value of A is this case [the ratio of cosmic variance to measurement error] is equal 
to 928/0.0124 — 75, 000. The Fisher one-sigma error on the fractional amplitude is a/2/5(1-|-1/A), essentially equal to 
the 0.71 cosmic variance limit. This large error hides the extreme non-gaussianity of the likelihood function. Indeed, 
even using I = 2 only, the statistical probability that there is no signal is infinitesimally small. To quantify this, 
consider the ratio of the likelihoods for two different models: (i) the best fit ACDM model with 2C2/7r ~ 928{^K)^ 
and (ii) a model where C2 — with no signal. Since the measured value of 3C2/tt is 201{fiK)^, the ratio of these two 
likelihoods - using only I — 2 data - is 



Ci _ 1 r_5 

£2^(1 + A)5/2 ^""P \ 2 



201/0.0124 201 



1 + A 



0.0124 



(22) 



of order 6"*"^°"°! Therefore, while the measured I — 2 moments give very little information about the amplitude of the 
anisotropies, they weigh in very heavily on the question of whether or not anisotropies have been detected. 

Returning to B-mode detection from the CMB, instead of using the one-sigma Fisher error, we compute the ratio 
of the likelihood for detection C{A = 1) vs. the likelihood of no detection £{A = 0). This likelihood ratio will vary 
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depending on the data. Fig.[7]shows the distribution of the hkehhood ratio C{A = 1)/C{A = 0) for 1000 mock "skies" 
(here sky simply means the five I = 2 CMB moments) generated from a true model with A — 1. Even when the signal 
to noise is unity (Ap = 1) so that the 1-sigma Fisher error is AA/A = 1.26, the likelihood ratio in a given experiment 
could be very large, signaling a detection. Ten percent of the mocks produced likelihood ratios greater than 100; that 
is, one would have concluded with 99% certainty that there are B-modes. A 90% CL detection would have occurred 
41% of the time. If the signal to noise were larger, Ap = 10, then detection is virtually assured, with a 99% detection 
emerging from 95% of the runs. 
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FIG. 7: Distribution of likelihood ratios for gravitational waves (GW) and no GW for 1000 fake skies generated when the true 
model contains GW's (using I = 2 CMB data only). When the underlying model has Ap = 1, the Fisher estimated 1-sigma 
fractional error on the amplitude of GW's would be y^2/5 x (1 + \p)/\p = 1.26, while the Ap = 10 model produces a Fisher 
error of 0.7. Nonetheless, a statistically significant detection is possible in the first case and virtually assured in the second. 

This suggests that, while lensing will not be of much use in constraining the GW amplitude, it might help in firming 
up the evidence for a detection. For example, suppose a CMB experiment measured a non-zero B-mode, but one 
wanted to compute the likelihood that this was due to GW or to foregrounds. With a CMB experiment only, the 
likelihood ratio of these two "models" (signal due to GW or signal due to foregrounds) would be unity: there would 
no way to tell them apart. How much would a lensing experiment improve on this? To answer this, I generated 1000 
realizations of both lensing and CMB I = 2 moments with Ap = 10 and Ai = 1 and a given value of a. For each 
realization, I computed the likelihood ratio of two models: 

• Model 1: True Model Signal to Noise in CMB Ap = 10; in lensing A^ = 1 and the true value of a 

• Model 2: Null Model Foregrounds produce Ap = 10, soAp = 

The distribution of these likelihood ratios is shown in Fig. [S] for several values of a. When a = 0, (7,39)% of the 
realizations led to a (99,90)% detection; For a = 0.5 those percentages go up to (12,45), and for larger values of 
a = 0.9 up to (21,57). So when the question posed is one of detection, as opposed to parameter determination, 
non-zero cross-correlation is beneficial. 

A harder question is whether one could detect the cross-correlation. To approach this question, I set the true values 
of Ap = 100 and Ap = 1 and computed the likelihood ratio of two models with: 

• Model 1: a equal to its true value 

• Model 2: a = 

The likelihood ratio was rarely significant. A 90% detection is never obtained if a = 0.5 and only 11% of the time 
even if a = 0.95. So we are unlikely to detect the cross-correlation of the lensing and CMB B-modes. 
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FIG. 8: Likelihood ratios of GW vs. (foregrounds in CMB and no signal in lensing) for 1000 mock skies generated from a true 
model with GW. 



VI. CONCLUSIONS 



Primordial gravitational waves produce indirect effects on both the polarization of the CMB and the lensing of 
distant galaxies. These effects are correlated, with a size and sign which depends on the angular scale and on the 
redshifts of the background galaxies. The correlation may help sort out systematics if a B-mode detection is made 
in the CMB. It seems likely that all-sky lensing surveys will be carried out for other purposes, so the search for the 
B-mode in that arena will cost nothing. It might even be used to motivate such surveys should the CMB B-modes 
be detected. 

I thank Anthony Challinor, Rob Crittenden, Wayne Hu, Lam Hui, Matthew Klcban, Eiichiro Komatsu, Hiranya 
Peiris, and Albert Stebbins for helpful conversations. This work was supported by the DOE at Fermilab and by NSF 
Grant AST-0908072. 
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Appendix A: Polarization Calculation 



Here I compute the transfer functions when the wavevector k lies along the z axis. Consider first the integral over 
n' in Eq. ([8]) for the case of a = + (the x contribution is identical so we need compute only this component). Thus, 



sin^ 9 cos(20) 



(r22(n')+ 5^2-2 («')) 



(Al) 



where 9 and 4> are the standard spherical angles with respect to the z-axis. Meanwhile the exponential can be expanded 
in spherical harmonics, or since k is taken along the z-axis: 



^^k■n'AD' ^ ^ ^ 471(21' + 1^' Yi,o{n')ji, (k AD') . 



(A2) 



Therefore, 



jW(H,Ai?') 



327r2 

"Is" 



l'=0 



J2 VWTTi''ji'ikAD') / d^n'Y*^„{h')Yi,o{n') {Y22{n') + Y2-2{n)) ■ (A3) 



The integral of the product of the Yim's over h' can be expressed in terms of the Wigner 3j-symbols 
(fn'Yl,„{n')Yi,o{n') {Y22{n') +Y2-2{n')) = 



25(2^ + 1) (2 2 1' 




An 



2 2 1' 
2 m" 



2 2/' 

-2 m" 



The first 3j-symbol is non-zero only when V = 0, 2, 4, and then is 



2 2 r 





\(V 



[{I'm 



112 



(4-0! 
(5 + /')! 



^'"^ {2 + 1' /2)\ 
{2~l'/2)V 



(A4) 



(A5) 



The first 3j-symbol in the square brackets is zero unless m" = — 2 while the second is non-zero only if m" = 2 and is 
equal to the first, so 



jit'(fcz,AD') = -A/^ 5^ i2l' + l)jAkAD')j- 



I'l 



i'=0,2,4 



[{I'm 



112 



(4-r)!]i/^(2 + r/2)! /2 2 I' 

(2-/'/2)!^""'"^ ""'^^^2 -2 



(5 + /')! 



(A6) 
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The relevant 3j-symbols with I' = 0,2,4 are l/\/5, ■y/2/35, and l/V 630 respectively. Therefore, 



407r 



35 



40^ 



+ (5m' 



m",2j 



ij2{kAD') 
[kAD'f ' 



(A7) 



A similar calculation for the integral in Eq. ^ begins by expanding the exponential in terms of spherical harmonics. 
We then are left with an integral over the product of terms of the form Yi'm.' (n)+?.Y-2m" (n)+2Y,t, (n). The first two of 
these can be re-expressed as a sum over ±2Ylm using Eq. (8) from Ref. [37j. The integral will pick out only the term 
with L = I and M = m, so in terms of the Clebsch- Gordon coefficients. 



Yl'rn'in)±2Y2rn"in) 

The integral is then 



5(2Z' + 1) 



_47r(2/ + 1) 



1/2 



(2, m", m'|2, 1'; I, m) (2, 1'; t2, 0|2, I'; I, ±2Yi,n{n) . 



Ii,nm"ikz, D) = %/5(5„,,„,» ^(2r + ji>{kD){2, 1'- to, 0|2, I'; I, 



V 2 I 
0-2 2 



V 2 I 
2-2 



(A8) 



(A9) 



where one of the Clebsch-Gordon sets has been written in terms of Wigner 3j-symbols and the identity y;,m-m"(^) — 
{-'iySm,m"V2l + I/a/Itt has been used. The second Wigner 3j-symbol is identical to the first apart from a sign flip 
in the bottom row, which changes the overall sign by (—1)' +2+'. Only if that sign is negative will cancellation be 
avoided, so one requirement is that I' + I is odd. The values of Z, I' for which these 3j-symbols are non-zero are those 
which satisfy: \l ~ l'\ < 2 < \l + l'\. The only values of I' which satisfy both conditions are I' = I ±1. Therefore, 



Il„irn"{kz,D) = 2\/5(5„ 



{21 - l){-iy-^ji^i{kD){2, l-l;m, 0|2, m) 



l-l 2 I 
-2 2 



+(2/ + l){-iy+'ji+i{kD){2, / + 1; TO, 0|2, l + m) 



l + l 2 I 
-2 2 



(AlO) 



Since we will be multiplying this by J^, which above was shown (in the case k \\ z) to vanish unless to" = ±2, we 
need evaluate the Clebsch-Gordon coefficients only for m = m" = ±2. For example. 



(2,?-l;2,0|2,/-l;Z,2) 



l-l 2 I 
-2 2 



(-i)'V2lTT 
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( iV+i /;Tr^ (2^-3)!(/ + 2)! 
^-'^ ^^\2l + 2)\{l-2)\ 



so the coefficient of for to = 2 is 



8V5(2. - l)H)-(-l)-V2m|^||±| = .-/Ilia + 2). 



(All) 



(A12) 



Using the symmetries of the 3j-symbols, it is easy to see that the to = —2 term is equal and opposite. Carrying out 
similar steps for the coefficient of leads finally to 



,{kz,D) 



21 + 1 



^m.m" \ym,2 



<5™,_2) [{I + 2)ji-i(kD) -{I- l)ji+i{kD)] . 



(A13) 



Therefore the transfer function in Eq. ([7]) reduces to Eq. ((T 



Appendix B: Weak lensing Calculation 



In this case, contracting the indices and expanding the exponential leads to 



T,t;(+)(fcz) = k dDT{k,rjo -D)f2 ^' ]J ^^^-^^^jAkD) ( <fnYC^{h)Yi,,{h)[Y22{h)-Y2-2{n)]. (Bl) 
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The integral over the product of spherical harmonics can again be expressed in terms of the Wigner 3j— symbols, 
leading to 



ni,( + ) 
'-Im 



{kz) 



2(2/ + l)7r 
3 

I V 2 
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(B2) 



Only I' = 1,1 ±2 are non-zero here, so the integrand will have terms proportional to ji,ji±2- The coefficients of these 
terms, however, are such that the resulting sum is proportional to ji{kD)/ (kD)^ , so 



Tt^'-\kz) = {^y-'[^m^^,-5. 
This leads directly to Eq. (fTTj) . 
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